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ABSTRACT 



Analysis of the hydrodynamic and helioseismic effects in the photosphere 
during the solar flare of July 23, 2002, observed by Michelson Doppler Imager 
(MDI) on SOHO, and high-energy images from RHESSI shows that these ef- 
fects are closely associated with sources of the hard X-ray emission, and that 
there are no such effects in the centroid region of the flare gamma-ray emission. 
These results demonstrate that contrary to expectations the hydrodynamic and 
helioseismic responses ("sunquakes") are more likely to be caused by accelerated 
electrons than by high-energy protons. A series of multiple impulses of high- 
energy electrons forms a hydrodynamic source moving in the photosphere with 
a supersonic speed. The moving source plays a critical role in the formation of 
the anisotropic wave front of sunquakes. 

Subject headings: Sun: flares - Sun: X-rays, gamma-rays - Sun: oscillations 



1. Introduction 



"Sunquakes", the helioseismic response to solar flares, are caused by strong localized 
hydrodynamic impacts in the photosphere during the flare impulsive phase. The helioseis- 
mic waves are observed directly as expanding circular-shaped ripples in SOHO/MDI Dopp- 
lergrams, which can be detected in Dopplergram movies and as a chara cteristic ridge in 



time-distance diagrams, (IKosovichev fc Zharkova 



calculating integrated acoustic emission (IDonea et al. 



1998 



. .Kosovichevll2006ar). or indire ctly by 
I999I : IPonea fc Lindsevl[2005[ ). Solar 



flares are sources of high-temperature plasma and strong hydrodynamic motions in the so- 
lar atmosphere. Perhaps, in all flares such perturbations generate acoustic waves traveling 
through the interior. However, only in some flares the impact is sufficiently localized and 
strong to produce the seismic waves with the amplitude above the convection noise level. It 
has been established in the initial July 9, 1996, flare observations (IKosovichev fc Zharkova 
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19981 ) that the hydro dynamic impact follows the hard X-ray flux impulse, and hence, the 
impact of high-energy electrons. Nevertheless, a common paradigm is that the sunquake 
events are caused by accelerated protons because protons carry more momentum and pen- 
etrate deeper into the solar atmosphere than electrons, which loose most of their energy in 
the upper chromosphere. This paradigm is not easy to test because the gamma-ray emission, 
which indicates the presence of high-energy protons, is rarely observed. 

In a large X17 flare of October 28, 2003, the gamma-ray emission observed by RHESSI 
was located close to the hard X-ray sources and two of the three places of the phototospheric 
impacts (sunquake sources) (IKosovichevI l2006al ) . Because of the close locations of the hard 
X-ray and gamma-ray sources these obse rvations could not exclude the possibility of the 
proton or mixed electron-proton impacts (jZharkova fc Zharkovl 120071 ) . 



However, in one event, X4.8 flare of July 23, 2002, the hard X-ray and gamma-ray 
sources were significantly separated from each other. The centroid of the 7-ray 2.233 MeV 
neutron-capture emission was found to be displaced by 20" ± 6" (with 5-sigma confidence) 
from that of the 0.3 — 0.5 MeV X-ray emission implying a difference in acceler ation and/or 
propagation between the accelerated electrons and ions (IHurford et al.ll2003l ). Therefore, 
this flare provides a unique opportunity to investigate the photospheric and helioseismic re- 
sponses separately for high-energy electrons and protons. In this Letter, I present results of 
the analysis of the relationship between the hard X-ray and gamma-ray emission s and the 



hydro dynamic and seismic signals in the photo sphere, using data from RHESSI (ILin et al. 



2OO2I ) and MDI on SOHO (IScherrer et al.lll995l ). RHESSI provides X-ray/gamma-ray imag- 
ing spectroscopy from 3 keV to 17 MeV with angular resolution 2.3" — 3' (35" at gamma-ray 
energies) over the full Sun. MDI measures the Doppler velocity and the line-of-sight mag- 
netic field of the photospheric plasma every minute with 2 arcsec/pixel resolution also over 
the full Sun. 



2. Analysis of SOHO/MDI and RHESSI data 



X4. 7 flare of July 23, 2002, was the first flare, for which gamma-ray images were obtained 
jHurford et al.ll2003h. _3ther examples of flares with RHESSI gamma-ray images are given 
by Hurford et al. ~ (j2006 ). The properties of the gamma-ray and hard X-ray emissions, and 
also other aspects of the July 23, 2002, flare, are discussed in the RHESSI special issue of 
ApJ Letters (v.595,no.2, 2003). The RHESSI observations revealed three hard X-ray sources 
and a gamma-ray so urce. Their positions in the MDI magnetogram are shown in Figure [1^ 
(iKrucker et al.ll2003l ). The hard X-ray sources (marked as fl, f2 and f3) were located on 
both sides of the magnetic neutral line. The morphology of the gamma-ray emission was 
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not resolved but it could not have been more than 1 arcminute (FWHM) in extent and its 
centroid was 20" ± 6 south from the centroid of the hard X-ray sources and about 30" from 
the fl source. 

The MDI Dopplergrams show strong impulsive variations close to the hard X-ray sources, 
but no impulsive variations in the region of the gamma-ray source centroid, or anywhere out- 
side the hard X-ray sources. Figure lb shows the positions of the impulsive Doppler velocity 
signals in the photosphere during the flare impulsive phase between 00:27 UT and 00:36 UT, 
July 23, 2002. This flare was rather close to the limb (coordinates of the flare sources are 
given in Fig. [1]); the distance from the disk center was approximately 70 degrees. Thus, it is 
possible that the projection effect contributed to the opposite sign of the Doppler shift across 
the neutral line, if the angle between the magnetic field lines directing the plasma motion at 
the footpoints (Fig. [2]) and the line of sight changes the sign. However, there might be other 
reasons related to flare hydrodynamics, which should be explored. 

The strongest Doppler signal, corresponding to a downward plasma motion, appeared 
near the X-ray source, fl. Its position moved during the impulsive phase in the North 
direction. This motion is discussed in more detail in the next section. The time dependence 
of the velocity signal at source fl corresponds very well (with the correlation coefficient 
of 0.8) to the total hard X-ray flux in the 50-300 keV range (Fig. |3t,d). The gamma-ray 
emission (Fig. [3^) is delayed by ~ 100 sec because of the time for the neutrons to thermalize 



flMurphv et al.l 120031 ) . 



The helioseismic waves are best visible at frequencies of about 5-6 mHz. To search 
for these waves the Dopplergrams were remapped into the heliographic coordinates, tracked 
to remove the displacement caused by the solar rotation, and then filtered using a band- 
pass filter centered at 5.5 mHz with a FWHM of 2 mHz. Then, the filtered Dopplergrams 
were remapped into the polar coordinates, centered at various points including all hard 
X-ray sources and the gamma source centroid region, and averaged azimuthally in several 
angular sectors. The averaged signals are plotted as a function of the radial distance and 
time, constituting time-distance propagation diagrams. The diagrams were inspected for an 
elongated characteris tic ridge-like structure, which is caused by helioseismic waves as pre- 
dicted by the theorv jKosovichev fc Zharkova 1995 ). and observed in other sunquake events 



( Kosovichev fc Zharkova 19981 : Kosovichev 2006a ). In this rather weak ridge ap 



peared only in the propagation diagram, which was centered in the region of the strongest 
impulsive Doppler signal at the fl source and averaged in the North- West quadrant. It can 
be identified in Fig. Ic for distances between 20 and 40 Mm, just above the theoretical time- 
distance relation for helioseismic acoustic waves. The traveling wave front can be seen in the 
movie of the frequency-filtered Dopplergrams. The observed signal is rather weak because a 
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flare gen erates high-frequency acoustic waves, in which the plasma velocity is predominantly 



vertical (iKosovichev &: Zharkovalll995l ). and its line-of-sight projection is reduced by almost 
2/3 due to the close-to-limb location. The amplitude of the line-of-sight plasma velocity in 
this wave was about 20 m/s. For the other central positions (including X-ray sources f2 and 
f3, and X-ray and gamma-ray centroids) and sectors, the seismic waves were not detected. 
This analysis puts the source of the seismic wave within the lower red part of the Doppler 
source fl in Fig. [T)d. The start time estimated from the theoretical time-distance relation in 
the ray approximation is 00:28-00:30 UT. It is interesting that X-ray source f2 is marginally 
stronger than source fl, but the impulsive Doppler signal is stronger in the fl position. 

The close correlation of the hydrodynamic and helioseismic responses with the hard 
X-ray flux source and the absence of any significant photospheric signal in the region of the 
gamma-ray centroid provide evidence that the source of the helioseismic waves is associated 
with the high-energy electrons and not with the high-energy protons. We note that while 
the RHESSI data do not exclude the presence of protons in the footpoints of the hard X-ray 
sources, for this conclusion it is important that the proton flux in the gamma-ray centroid 
area was not weaker than that in the HRX sources. This assumption is supported by the 
RHESSI data. For further studies, it would be important to put precise hmits on the proton 
flux at the X-ray footpoints and estimate the relative energetics of proton and electrons from 
RHESSI data. 



3. Moving hard X-ray and sunquake sources 



A characteristic feature of the seismic response in this flare and several others (IKosovichev 



2006al jb[) is anisotropy of the wave front: the observed wave amplitude is much stronger in 



one direction than in the others. In particular, the seismic waves excited during the October 
28, 2003, flare had the greatest amplitude in the direction of the expanding flare ribbons. 
The wave anisotropy was attributed to the moving sou rce of the hydrodynamic impact, 
which is located in the flare ribbons ( Kosovichev 2006a 3). The motion of flare ribbons is 



often interpreted as a result of the magnetic reconnection processes in the corona. When 
the reconnection region moves up it involves higher magnetic loops, the footpoints of which 
are further apart. Of course, there might be other reasons for the anisotropy of the wave 
front, such as inhomogeneities in temperature, magnetic fleld, and plasma flows. However, 
the source motion seems to be quite important. 

It is interesting that in the case of the July 23, 2002, flare the seismic source identifled 
in MDI Dopplergrams as a place of strong Doppler shift in region fl was moving mostly 
along the flare ribbon, and consequently the seismic wave had the strongest amplitude in 



- 5 - 



the direction close to the direction of the source motion (but not precisely; in this case, in 
addition to the other factors, stronger foreshortening on the East side might have contributed 
to the signal loss in the NE quarter). The Dop pler source mot i on nic ely corresponds to the 
motion of the hard X-ray source discovered by iKrucker et al.l (120031 ) . Figure shows the 
evolution of the hard X-ray sources, fl, f2 and f3, positions o n the magnetogram; and Fig.|DD 
shows propagation diagrams for this sources determined by lKrucker et al.l (120031 ). From the 
top panel of Fig. |3)d, the hard X-ray source, fl, traveled approximately 7 Mm in 5 min; 
this corresp onds to the mean sp eed of approximately 20-25 km/s. The maximum speed 
according to iKrucker et al.l (120031 ) reached 50 km/s. 



Using the MDI Dopplergrams, a similar time-distance propagation diagram was con- 
structed for the plasma photospheric velocity along the line of motion of source fl. Figure 
Hb shows the Doppler velocity along a 2-pixel wide strip along this line. This diagram shows 
that the evolution of the hydrodynamic impact source is very similar to the evolution of the 
hard X-ray source (top panel in Fig. Ht). The mean speed of the hydrodynamic source was 
also about 20-25 km/s. 

Therefore, we conclude that the seismic wave was generated no t by a single impulse 



as wa s suggested in th e sunquakes models of iKosovichev &: Zharkoval (119951 ): iMedrek et al. 



(I2OOOI ): iPodestal (120051 ) but by a series of impulses, which produce the hydrodynamic source 
moving on the solar surface with a supersonic speed. The seismic effect of the moving 
source can be easily calculated by convolving the wave Green's function (the wave signal 
from a point 5-function type source), G{x — Xs,y — ys,^) with a moving source function, 
S{xs — Vxt, ys — Vyt, t). The results of these calculations are illustrated in Fig. [5l which shows 
the wave front for a source moving along the x-axis with a speed of 25 km/s. The strength 
of this source varied with time as a Gaussian with FWHM of 3 min (it is shown by black 
diamond s). The Green's function was calculated by using the standard mode summation 
method (IKosovichev fc Zharkovalll995l ). The strong anisotropy of the seismic wave is evident. 
Curiously, this effect is quite similar to the aniso tropy of seismic wave s on Earth, when the 



earthquake rupture moves along the fault (e.g. iBen-MenahemI Il962l ). Thus, taking into 
account the effects of multiple impulses of accelerated electrons and moving source is very 
important for sunquake theories. These effects will be discussed in more detail in our future 
publications. 



4. Discussion 



The analysis of RHESSI X-ray and gamma-ray images and SOHO/MDI Dopplergrams 
of the July 23, 2002, X4.8 solar flare revealed that the hydrodynamic and seismic responses 
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are closely associated the hard X-ray emission, both spatially and temporally, but showed 
no significant responses in the gamma-source centroid area. Because this fiare was one of 
strongest gamma-flares, and the hard X-ray and gamma-ray sources were separated, these 
observations show that the accelerated protons are unlikely to be a source of the hydro- 
dynamic response and sunquakes. Furthermore, th e detailed a n alysis of the dynamics of 
sunquake sources in this Letter and in the paper by iKosovichevI (j2006al ) reveals their close 
association with expanding flare ribbons and rapid HXR source motion along the ribbons, 
and, thus, with the magnetic reconnection process. The fast motion of these sources results 
in strong anisotropy of the seismic waves, clearly observed in the MDI data. 

The general picture that comes from the analysis of MDI and RHESSI data is consis- 
tent with the previously developed hydrodynamic thick-target model, illustrated in Fig. [2] 
jKostiuk fc Pikelnerl[l975l : llivshits et al.lll98ll : iFisher et allllQS^ lKosovichevl[l986l ). In this 
model, high-energy electrons heat the upper chromosphere to high temperatures generating 
a high-pressure region, expansion of which causes evaporation of the chromospheric plasma 
and a high compression shock. The shock reaches the photosphere and excites the seismic 
waves. However, the new results show that it is important to include effects of the multiple 
impact and moving source in the thick-target and sunquake models. 

The photospheric and helioseismic effects observed during the impulsive phase of solar 
flares are closely related to the processes of acceleration and propagation of electrons and 
ions, and may provide new important information about these processes. 
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X4.8 Flare, 23 July 2002 Flare (0-90 deg) 
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Fig. 1. — a) RHESSI observations of HXR (fl, f2, f3) sources (blue contours show the 50-100 
keV map with 3" resolution, and the blue cross shows the HXR centroid) and a gamma-ray 
source (orange circle shows the location of the gamma-ray centroid with la error) (from 



Krucker et al.l (120031 )). The gray-scale background is MDI magnetogram (±600 G range). 



b) MDI observations of Doppler velocity sources. Orange contour lines show positive (red- 
shift) velocity greater than 0.6 km/s; the red contours show 1 km/s; the blue contour lines 
show negative (blue-shift) velocity of -0.6 km/s, and dark blue shows -0.7 km/s. c) Time- 
distance map revealing a seismic wave front, which travels in the North- West direction from 
the location of source fl. The yellow dashed curve is a theoretical time-distance relation for 
helioseismic acoustic waves. 
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Fig. 2. — Illustration of the hydrodynamic thick-target model and the mechanism of sun- 
quakes. High-energy electrons accelerated in the upper corona are injected along magnetic 
field fines into the atmosphere, generate a hard X-ray emission in the loop footpoints and 
heat the upper chromosphere to high temperature, producing a high-pressure region. The 
high-pressure region expands producing upward and downward propagating shocks. The 
downward shock reaches the photosphere and causes a sunquake. 
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X4.8 flare, July 23, 2002 
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Fig. 3.- 
of time. 



Integrated X-ray and gamma-ray fluxes, and Doppler shift in source fl as a function 
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Fig. 4. — a) Evolution of HRX sources. The increasing size of the symbols represents times 
from 00:26:35 to 00:39:07 UT, and b) HXR profiles (black in the top two panels shows 
enhanced emission) along the ribbons, showing motion with speed of up to 50 km/s (from 



Krucker et al.l (120031 )). c) Doppler velocity profiles along the fl ribbon, showing motion with 



averaged speed 25 km/s. The inclined orange lines correspond to 25 km/s (long dash) and 
50 km/s (short dash). 
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Fig. 5. — Theoretical model of seismic waves from a moving source, which explains the 
observed anisotropy of sunquakes. The point impulsive source is moving in the x direction 
with the constant speed of 25 km/s. Its strength as a function of time has a Gaussian shape 
with FWHM of 3 min. The locations of the source are shown by black diamonds at the 
center, the size of which is proportional to the source strength. 



